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The energy arising from the hyperfine interaction, the interaction with the 
magnetic field and the cross term between them have been studied for ortho H 2 
at c 3 7 t u state. The cross term gives rise to a first order correction of the g F 
values which accounts for the discrepancies between the theoretical zero th order 
g F and the experimental values. The interactions are expressed as products of 
irreducible spherical tensors and Racah's method is used to derive all the diag- 
onal and off diagonal matrix elements. The case of rotational level, N = 1 , is 
presented in detail. The coupling constants of hyperfine structure have been 
evaluated by using a single configuration electronic wave function for c 3 77 u 
state of H 2 . //-fojhcte 


INTRODUCTION 


The odd rotational levels of ortho H 2 and the even rotational levels of para 
H. of the c 3 77 state are metastable with a life time of the order of milliseconds. 

2 u 

Licnten 1,2 measured the fine structure splittings and the g^ values of the rota- 
tional level N = 2 of para H 2 to a high precision by atomic beams magnetic 
resonance method. The hyperfine splittings and the g values of the N = 1 level 
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of ortho H 2 have recently been measured by Lichten and Brooks. 3 Since hydro- 
gen molecule is the lightest diatomic molecule, pure Hund’s case 4 b coupling of 
angular momenta can be applied to a very good approximation. The properties 
of the c 3 7t u state has therefore been of great theoretical interest. The fine 
structure of this state has been studied by Fontana 5 and by Chiu. 6 Its hyperfine 
structure has been studied by Frey and Mizuschima, 7 and the relative life time 
of the fine structure levels due to the perturbation of other electronic states 
has been studied by Chiu. 8 

The fine structure separations are of the order of 5000 Mc/Sec 2,6 which 
is not very large in comparison with the hyperfine separations (of the order of 
500 Mc/Sec). The second order corrections to the hyperfine energies will 
therefore be large. Recent measurements 3 indicate that the measured g F values 
differ considerably from the zero th order (g£ 0) ) values where the perturbations 
from the nearby fine structure levels have been neglected. These differences 
are mainly due to a nonvanishing cross term of the hyperfine interaction and the 
interaction of the molecule with the magnetic field. The corrections to the g F 
values from this cross term will be of the first order (g^ 1 ) ) , and they can take 
on large values. 

In section 1 of this paper we will discuss these corrections on g F values. 

The diagonal and off diagonal matrix elements of the hyperfine interacting 
Hamiltonian will be derived in section 2. In section 3 we will work out the 
details for the case of rotational level N = 1, namely the second order energy 
corrections E (2) and the corrections on gp will be given for all the hyperfine 
states. A set of coupling constants which have been evaluated numerically from 
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a single configuration electronic wave function 6 for the c 3 ^ state will be given 
in the last section. 

1. ENERGY CORRECTIONS AND g p CORRECTIONS 

The perturbing Hamiltonian, H', of a diatomic molecule due to its hyperfine 
interaction (H^ s ) and its interaction with the external magnetic field (H^,) is 
the following: 


H' = H hfs + H m , 


where 


H hfs = Hj (hfs) + H 2 (hfs) + H 3 (hfs) 


and 


E 1 (hfs) 


(..-•/2)r (I -.h) 

cj 1 


i 


r 


3 

i/3 


H 2 (hfs) - (477g t a a /3)/ t S(r i/3 ) S i • Ig 


H 3 (hfs) - - (g r a 2 / 2 


fe • .y - 




r . 


i/3 


J r i/3 


( 1 - 1 ) 


The indices i ( = 1, 2) refer to the two electrons and the indices / 3 ( = a, b) refer 
to the two nuclei. and Sj are the orbital angular momentum and the spin of the 
electron i ; I^ is the spin of the nucleus / 3 ; g t = 0.003038 for proton and a is the 
fine structure constant. 
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H w is the interaction of the magnetic moments of the molecule with the ex- 
ternal magnetic field, 

H m = MoUsS + g^N + gjl) • H (1-2) 

where /x Q is the Bohr magneton; S is the total electronic spin, S = Sj + S 2 ; N 
is the total rotational angular momentum which is the resultant of pure rotational 
angular momentum 0 and the projection of electronic orbital angular momentum 
along the figure axis A; I is the total nuclear spin, 1=1, + I„ ; and H is the 
external magnetic field. J5 and N will couple to give the total electronic angular 
momentum J , J = N + S , and J will couple with I to give the total angular mo- 

v*'"'- VvV l - v ' 1 ' vV ' / . "i — V V-vA- 

mentum JF of the molecule, F = J + I . 

Due to the closeness of the fine structure separation, there are mixings 
among the levels of same JF but different J’s. In a zero magnetic field, the 
wave function of a hyperfine level (F, m F ) within a fine structure level J (for a 
given rotational level N of c 3 7 t u state) is expressed up to first order correction 
as, 




J.F,m F 


V (X’FJiNlHhJ VFJ'IN) 

0 JF,m F + L_^ Ej p - E . F ^J'F,m F (1-3) 

J' 


where <£j jF is the ortho normal zero th order wave function for level QF, m F ) 
of the c 3 7t u state where both F and J are good quantum numbers. In the pres- 
ence of a very weak magnetic field the energy of this hyperfine level is therefore 
(up to second order correction) the following, 
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E 


J.F.m F 


Ej.F ( hf S ) + (g F (0> + %F ) ^O m F H 


(1-4) 


where 


Ej,f ( hfs ) = ( 37T u ,FJIN|H hfs | ^FJIN) 



( 3 7T u ,FJIN|H hf J 3 7 t u ,FJ'IN) 


J'*J 


E j,f E j',f 


(1-5) 


The first and the second terms on the right hand side of Eq. (1-5) are the first 
order energy, E^ (hfs), and the second order energy, E ( j 2 ^ (hfs) , respectively 
due to the hyperfine interaction and, 

g F (0) (J) = [F(F + 1)]- 1 /2(FJINS II g N N + g s S + gj. Ijl FJINS) (1 _ 6) 


g F ' (J) = 


2 

[F(F + 1)] 1/2 



J'^J 


( 3 77 u ,FJIN|H hfs | 3 7 t u ,FJ , IN)(fJ' INS || g N jj+g s S 1| FJINS) 

E J,F _ E j.'f 



1 1 (X.FJINiHhJ 3 v u ,rj'lN) 


J't^J 


E J,f E J',f 


g F (0 ) (J') (1 ' 7) 
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g< F °) (J) is the zero th order g p value when both F and J are good quantum 
numbers. The first term on the right hand side of (1-7) is the first order cor- 
rection, g^ 1 ) (j) , which can be of comparable magnitude to that of g ( F 0) (J) 
due to the energy separation in the denominator is not large. The second term 
in the RHS of Eq. (1-7) is the second order correction, g(. 2 > (J) , which is small. 
The reduced matrix elements in Eq. (1-6) and Eq. (1-7) are defined by Wigner- 
Eckart'theorem 10 as follows, 

(F', m' |T/ l >| F, m) = C(FLF; m/jm ' ) (f' || T< L > II F) 

where T a) is the component of the L th rank irreducible tensor T ( ' L) , and 
C(FLF'; m/um') is the Clebsch-Gordan's coefficient . 11 

2. MATRIX ELEMENTS 

To obtain the hyperfine energy separations and the corrections on g p values, 
one needs to evaluate the diagonal and the off diagonal matrix elements of the 
interacting Hamiltonian H'. For Hand's case b diatomic molecule, the angular 
momentum N, S and I are well defined with respect to the space fixed coordinate 
system. The matrix elements of H„ over the molecular wave function can there- 
fore be easily evaluated. 12 H x (hfs) is the sum of the scalar products of the 
electronic angular momentum and the nuclear spin I For matrix elements 
diagonal in total nuclear spin I, one can write H x (hf s) in the following form, 

»i (hfs) = ( gl a 2 /4) (*' > ‘ JL 
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where x' refers to the fixed space coordinate system. By considering the above 
expression as a scalar product of two tensors of first rank, one can evaluate the 
matrix elements of Hj (hf s) as follows, 



, FJ' IN 



*&. (x' ) • I 

VWV A V WV 7 vW 


3 vr u , F J I N 


) 


= (-1) I+ J'- F [( 2J ' + 1) (21 + 1)] 1/2 W(J I J' I; FI) 


X 




where W(J I J' I; FI) is the Racah’s coefficient, 13 and 1=1 for ortho H 2 . The 
reduced matrix element (I | [l| |l) = [1(1 + 1)] 1/2 = /2 . After (1/ r?^) •£. (x/) 
been transformed from fixed coordinate system (x/ ) into molecular moving 
coordinate system (x) , the reduced matrix element of (l/r 3 ^ ) ^ (^J) is evalu- 
ated as follows. 14 
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where <((l/ r f z )> is to be evaluated in the molecular system (x) over the 

electronic wave function of c 3 7t u states, the axis z is the figure axis. The matrix 
element of Hj (hf s ) is therefore, 

( 3 tt u , FJ'IN |Hj (hf s ) | 3 t t u> FJIN ) 

= ( - 1) N+1 ~ F [6( 2J ' + 1) ( 2J + 1) (2N+ 1)/ N(N + 1)] 1/2 

x W(J 1 J ' 1 ; F 1) W(NJNJ' ; ll)a, 

where j' = J , J ± 1 and a is the coupling constant in terms of the integral 
<^(l/r 3 ^ ) l. iz y and it is defined in Eq. (2-9) below. 

H 2 (hf s) is the contact interaction of electron spin with nuclear spin. It can 
be written into an equivalent form for matrix elements diagonal in I and S , 

H 2 (hfs) = (^gj a2 /3) ^ S ( r i^)? ‘ .1 

i./S 


The matrix element of H 2 (hf s) is similarly decomposed into a product of Racah's 
coefficient; and reduced matrix elements, 

(F J' I N S | S • I | F J I N S) 

_ (.^N-s-F j-^j' + !) (2J + 1) (21 + 1) ( 2S + 1)] 1/2 

x W(JIJ ' I; FI) W(SJSJ ' ; Nl) (i || I || i) (s || S.|| s) 
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with 1 = 1, S = 1, (l|| l||l||) =/2 and (S||s||S) =/2 then, 

( 3 tt u , FJ' N |H 2 (hfs) | 3 tt u , FJN ) 

= (-1) N_F_1 6 [( 2J ' + 1 ) (2J + 1)] 1/2 W( J1J ' 1; FI) 

X W( 1J 1J ' ; Nl)b 

where b is in terms of the integral ^S(r i/8 )^ over the electronic wave function 
and it is defined in Eq. (2-9) below. 

H 3 (hfs) represents spin-spin interaction between electrons and nuclei. 
For matrix elements diagonal in S and I , it is equivalent to write H 3 (hf s) as, 



Since the doubly degenerate -n state is a linear combination of A = + l state and 
A = - 1 state, and the matrix element of H 3 (hfs) is nonvanishing between state 
of A = + 1 and state of A = - 1, we therefore have, 


( 3 77 u , FJ ' IN I H 3 (hf S ) I 3 7T u , FJIN) 


= (FJ' INS, A = 1 |H 3 (hfs) IFJINS, A = l) 

- (FJ'INSAn 1 | H 3 (hfs) IFJINS, A = - 1), 


(2-1) 


for odd rotational level (N = odd) of ortho H 2 - 15 If one expresses H 3 (hf s) as a 
sum of H 3 1/3 (hfs), 


H 3 (hfs) 


J2 ^ < hfs > 

i.P 
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where (hf s) is a scalar product of two first rank tensors and its matrix 
element can again be decomposed into a product of Racah's coefficient and 
reduced matrix elements. 

(FJ'INSA' |H 3 i/ 3 (hf S )| F J I N S A ) 


= _ (gj a 2 /8 ) (-1) i+J '- F [(2I + 1) (2J'+ 1 )] 1/2 W(IJIJ';Fl)(l || I || i) 


J ' N S A' || (7V) [S - (S * ] HJNSA 


where (I || 1 1| I) = [1(1+ 1 )] 1/2 and 1 = 1 for ortho H r Let, 


r (l) 


l- 3 ii/j (§.•+«>)/' 


i/3 


where x n > is a spherical tensor of the first rank and it can be constructed from 
a first rank spherical tensor and a second rank spherical tensor c£ 2 ^ iy g) 

where = (47t/5) 1/2 (£'9') • We therefore have , 10 

x m ( = k C( 1215 m ~ ^ m) S ^ (1> C -»-? 

A* 

where C( 121; /i.nu^m) is the Clebsch-Gordamcoefficient , 12 and 9*^ and 9 '.^ 
are the polar and azimuth angles of j- i(S with respect to space fixed coordinate 
system. The constant k = /10 can be obtained by considering the case of m = 0 . 

The reduced matrix element of x< are then decomposed into the product of 
reduced matrix element of C (2 > (0^ 9 ' ) and of S, and a 9j symbol. By using 
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the relation of (S||sJ|S) = [S(S + 1)] V2 and S = 1 for triplet electronic state, the 
reduced matrix element of may be evaluated as follows, 16 


(J'N S A' HxJ^/rJg || JNSA) 


= [l80 (2J + 1) ( 2N + 1)] 1/2 



1 

N 

J 


2 t^NA' || ^ )( - e i{3Vi{3> 



(2-3) 


The above 9j symbol is a sum of the products of Racah’s coefficients. 163 ’ 17 For 
the case of J' = J, J ±1 the 9j symbol can be expressed into the following sim- 
ple expressions, 



[4 +N(N + 1) ~ J(J + 1)] [3J(J + 1) +N(N + 1)] - 12 
12[SN(N + 1) (2N+ 1) (2N + 3) (2N- 1) J(J + 1) (2J + 1)] 1/2 



1 

2 

1 


[(J + N + 3) (N + J) (J - N + 2) (N - J + 1)] 1/2 [3( J + l) 2 - N(N + 1) - 6] 
12 [5N(N + 1) (2N+ 1) ( 2N + 3) ( 2N - 1) (J + 1) (2J + 1) (2J + 3)] 1/2 


After a proper transformation of the fixed coordinate (O' <p' n ) into the moving 

v i/3 T i/3 y - 

molecular coordinate ,9.^) we have, 14 

(nA' l|C< 2 >(^ /S cp' /S )/rf /S || NA) = (-l) k C(N2N; A'kA)(A'|c£ 2 )(^9 i/Q )/rf /S |A) 

(2-5) 
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where (<9.^ cp.^) /rf^y is independent of the molecular rotation and is to be 

evaluated over the electronic wave function of the c 3 7t u state of H 2 . By substi- 
tuting (2-2), (2-3), (2-4) and (2-5) into (2-1) one obtains the matrix element for 
H 3 (h f s) in the following, 

( 3 *V FJ'N|H 3 (hfs)| 3 77 u , FJN) = 


(-1)J'+ F W(1J1J'; FI) 


270 N(N + 1)(2N + 1) (2 J * + 1)(2J + 1) 
(2N - 1) (2N + 3) 


1/2 


3 - N(N + 1) _ /6 

N(N + 1) 2 


r \ 1 
| N N 

J' J 



(2-6) 


where J' = j, j ± 1, and the coupling constants c and d are defined in (2-9). 
Combining the matrix elements for H } (h f s ), H 2 (h f s) and H 3 (h f s) we have the 
matrix elements for H hfs as follows, 


K.F, J,N|H hfs | 3 tt u ,F, J,N) 


r 

Qi 1 fj(j + l)'+ N(N + 1) - 2] 


2J (J + 1) 


N(N + 1) 


f4 + N(N + 1 ) ~ J(J + 1)1 [3J(J + 1)+N(N+ 1)] - 12 
+ [2 + J(J + 1) - N(N + l)]b+ 2(2N - 1) (2N + 3) 


X 


N(N + 1) - 3 yg 
N(N + 1) c + 2 


d 


- 1 j 


(2-7) 
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( 3 ^ u ,F, J + l,N|H hfi | 3 vt u ,F, J,n) 


Q 2 Q 3 f a [3(J + 1) 2 - N(N + 1) - 6] 

2(J + 1) |N(N + 1) + b + 2(2N- 1) (2N + 3) 


x 


N(N + 1) - 3 
N(N + 1) c 



(2-8a) 


and 


(X,F, J - 1, N |H hfs | 3 7 t u ,F, J,N) 

Q 2 Q 3 f a [3J 2 - N(N+1) - 6] 

" " 2J |N(N + l) +b+ 2(2N- 1) (2N+3) 


N(N + 1) - 3 

N(N + 1) c + 2 d 


(2-8b) 


where a, b, c and d are the coupling constants to be evaluated over the electronic 
wave function of the c 3 ta u state and they are defined in the following: 

■ 2 4 , 


Bi 


a 


lz v 1 z 2 z 2z 

3 " + 7T + 7T + “ r 3 

la lb 2a 2b/ 


b - ( gl 7ra 2 /3) <S(rJ +8(r lb ) + S(rJ + § ( r 2b )) 


gj; ^ 


C 



\ 


la 


.) C< 2 > (0 lb <p lb ) C< 2 > (0 2a <p 2a ) C< 2 > (0 ab <p ab )_^ 


lb 


„3 
A 2a 


2b / 


d = 


«. «' /ci*> C<»(9 lb ?,„) C<*> Cl*> <p 2b ) \ 


la 


lb 


2a 


2b 
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(2-9) 



and 
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The coefficients Qj, Q 2 , Q 3 , Q 2 and Q 3 appeared in (2-7) and (2-8) are defined 
as follows, 

Q, = [F(F + 1) - J(J + 1) - 2]/2 

Q 2 = [(F+J+3) (F + l-J) (F + J) (J-F + 2)] 1/2 /2 


(J+N + 3) (J+N) (J-N + 2) (N-J + l) 
( 2 J + 1 ) ( 2 J + 3 ) 


1/2 


Q 2 - [(F+J + 2) (F + 2-J) (F + J-l) (J-F + l)] 1 ^ 


(J+N + 2) (J+N-l) (J-N + l) (N-J + 2) 


(2J-1) ( 2 J + 1 ) 


1/2 


(2-10) 


3. THE CASE OF N = 1 


With all the matrix elements developed in the last section, we will now ex- 
press the hyperfine energy and the g p values for the rotational level N = 1. 


S JiF (hfs) = E/y (hfs) + E/V (hfs) 


Substituting (2-7) and (2-8) into (1-5) we have, 


Ej<y (hfs) = (Qj/40) jlOa + 20b + [l6 - 3J(J + 1)] [V6d - c]J (3-1) 
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and 


1 


where 


E/.v (hfs) 


~ Q 2 Q3P] 2 
40(J + 1>J AE JJ+1 (F) 


+ 


Q 3 p 

40J 




AEj J±1 (F) = E(J, F) -E(J ±1, F) 

P S 10a + 20b + [3(J + l) 2 - 8] [V6d - c] 


(3-2) 


(3-3) 

(3-4) 


and 

P = 10a + 20b + [3J 2 - 8] [/6d - c] 

The Q’s and a, b, c, d are defined in (2-10) and (2-9). 

For gp values, we neglect the small correcting term g< F 2 ) (second term) 
in Eq. (1-7), then, g^. (J) = g(.D (j) and 

g F (J) = g K ( 0 ) (J) + g F (1) (J) 


For the case of N = 1 , we have 

g F ( 0 ) (J) = 


F(F+I) + J(J + 1) - 2 F(F + 1) - J(J + 1) + 2 

“ 4F(F + 1) (8s + 6h) + 2F(F + 1) «i 


g F ( (J) " 40F(F + 1 ) 


O 2 O 2 P 

'-Z 


0 2 0? P 

~ J - o 




/o n \ 
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where P and P are defined in Eq. (3-4), and g N = A/N(N+1)= 1/2, g s = 2 and 
gj = 0.003038.- Since the wave function *Pj FM is normalized up to the first order 
(see Eq. (1-3)), we therefore have 18 


2]gF (1) (J) = 0 . (3-7) 

J 

All the values of g(°) (J) for different N's have been tabulated 1 by using Eq. 
(3-5). We now summarize all the g ( F 1:) (J) and E< 2 ). (hfs) for the case of N = 1 
in the following, let 


K, 

k 2 

we have, 

g F (1) (F = 1, J = 0) 

i F (1) (F = 1, J = 1) 
g F (1) (F = 1, J = 2) 
g F (1) (F = 2, J = 1) 
g F (1) (F = 2, J = 2) 
g F (1) (F = 3, J = 2) 


= a + 2b - (/6d-c)/2 
= a + 2b + 2(/6d -c)/s 


5K 2 /16AE 12 (F = 2) 
-5K 2/ /l6AE 1>2 (F = 2) 
g F (1) (F = 0, J = 1) 


(3-8) 


(F = 1) 


0 (3-9) 


-5K 1 y / 3AE lj 0 (F = 1) 

25K 2 /48AE 12 (F= 1) + 5KJ/3AEJ0 
-25K 2 /48AE 12 (F = 1) 
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and 


E< 2 > (F = 1, 

J = 0) 

= = 1) 


E< 2 > (F = 1, 

J = 1) 

= 5(K 2 ) 2 /48AE 1j2 (F = 1) + (K x ) : 

2 /3AE lt0 (F = 1) 

E (2) (F - 1, 

CN 

II 

•— > 

= -5(K 2 ) 2 /48AE 1j2 (F = 1) 


E< 2 > (F = 2, 

J = 1) 

= 3(K 2 ) 2 /l6AE 1>2 (F = 2) 


E< 2 > (F = 2, 

J = 2) 

= -3(K 2 ) 2 / /l6AE 1 2 (F = 2) 


E< 2 > (F = 3, 

J = 2) 

= E< 2 > (F = 0, J = 1) = 0 

(3-10) 

Where the AEj j < (F) ’s 

are defined in (3-3). Note that the sum rule relation 


(3-7) holds for the g^ 1)! s in (3-9) and a similar relation E j 2 p = 0 also holds 

j 

for E (2) in (3-10). 

Experimentally, one can determine the constants and K 2 from the fine 
structure and the hyperfine structure separations and the g j, 1 ^ values. Knowing 
and K 2 one can obtain the linear combinations, (a + 2b) and (/6d - c) , from 
Eq. (3-8). Measurements on N = 1 rotational level alone therefore cannot deter- 
mine the coupling constants a, b, c, and d independently. Additional measure- 
ments on N = 3 level will be needed to determine these hyperfine coupling 
constants independently. 
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4. COUPLING CONSTANTS 


The coupling constants a, b, c, and d are defined in (2-8). They have been 
evaluated by using a single configurational electronic wave function 6 of the c 3 tt u 
state. The wave function is a linear combination of the Heitler-Condon type with 
coefficient C H and ionic-type wave function with coefficient C . Both types of 
wave function are made of Slater-type atomic orbitals of lscr and 2p77 symmetry. 
For the orbital exponents £ lsCr (= 1.24108) and £ 2p7T (= 0.50282) the coefficients 
Cj (= 0.2808) and C H (= 0.5612) are the values which minimize the total energy 
E ( = -0.718793 a.u.) of the molecule in a c 3 tt u state at equilibrium internuclear 
distance R e (= 1.96080 a.u.). The values of the coupling constants are evaluated 
numerically 19 on the IBM 7094 computer, and they are the following. 

a = 17.44 Mc/Sec, b = 545.27 Mc/Sec. 

c = 28.95 Mc/Sec, d = -8.27 Mc/Sec. 

Frey and Mizushima 7 have also computed these coupling constant by using 
Amemiya's wave function. 20 Their c and d. values appear to be very different 
from ours. 21 It seems that they may have defined their constants differently, 22 
although the definition of coupling constants was not given by them. 

With the calculated coupling constants one can obtain all the first order 
gU) » s from Eq. (3-9) and the second order E (2) 's from Eq. (3-10). The g F values 
thus corrected by g ( F * > *s give a better agreement with the experimental values 3 
than the uncorrected g^ 01 * 's. Since all the operators involved in calculating the 
coupling constants are one electron operators, a set of much improved coupling 
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constants can therefore be calculated easily when an accurate electronic wave 


function for c 3 7 t state is available. 

u 
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and it is incorrect. From their Eq. (5) and Eq. (6), 


"d " is /6/4 times of what we have defined here in Eq. (2-8), and this 
definition is inconsistent with their Eq. (10). Furthermore their Eq. (10) 
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appears to be incorrect because the ratio of the terms with c and d should 
be J independent whereas their ratio is J dependent (and this is the reason 
that c and d cannot be determined independently by measurements on 
rotational N = 1 level alone). 
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